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We consider the one-loop effective action due to gravitons in a FLRW background 
witli constant e = —H/H^. By expanding around e = (corresponding to an expan- 
sion around de Sitter space), we can study liow tlie deviation from de Sitter space 
effects tlie quantum corrected Friedmann equations. We find that, at zeroth order 
in e, one-loop effects induce only a finite shift in the coupling constants. At linear 
order in e there is however a divergent contribution to the equations of motion. This 
contribution leads to a nontrivial term in the renormalized equations that depends 
logarithmically on H and thus cannot be absorbed in local counterterms. We find 
that deviations due to this term are unobservably small. Our study shows that 
quantum effects in quasi de Sitter space can be fundamentally different then in de 
Sitter space, albeit in the case under consideration the effect is unobservably small. 

I. INTRODUCTION 

Because of the potential relevance for inflationary cosmology, the quantum behavior of 



gravitons on a (locally) de Sitter b ackg round has 

m,m B,flB BBQ HHQ lis 



years 
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3een a wide 

y,y [17 
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y,M y, 



21 



22 



23 



le past 



24|. 



One line of research deals with the back-reaction of gravitational waves on the background 
spacetime [l3|, [ij, [l7|, [l8|. However of more interest for the present work is the one loop 
back-reaction by virtual gravitons on a de Sitter background which has been calculated by 

25||26j. Since it is not clear whether in these 



several authors using different techniques [1 



works exactly the same quantity is calculated and the renormalization schemes differ, the 
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numerical coefficients differ. However the main result is qualitatively the same: one loop 
graviton contributions to the expectation value of the energy momentum tensor result in a 
finite, time independent shift of the effective cosmological constant. Since the contribution 



can always be absorbed in a counterterm 25|], the exact numerical coefficient has no real 

physical meaning. 

The goal of this paper is to go beyond the works mentioned above and calculate the one 

loop effective action induced by gravitons in a more general background space-time using 

dimensional regularization. The geometry we consider is a Friedmann-Lemaitre-Roberston- 

Walker (FLRW) geometry with Hubble parameter H = a/ a and the additional constraint 

that 

is a constant. Standard matter, radiation or dark energy dominated universes all satisfy this 
constraint (recall that in matter era e = 3/2, while in radiation era e = 2) and de Sitter 
space is the special limit when e -^ [27]. One immediate problem with working in such a 
space-time, instead of in de Siter space, is that, for consistency of the Einstein equations, 
the addition of matter fields is unavoidable. Whereas in de Sitter space, the only relevant 
metric fluctuations are the tensor modes (gravitational waves), in a more general setting also 



the scalar modes, due to 
to be taken into account 



;he mixing of gravitational and matter degrees of freedom, have 



28||29| 



21 



22 



23 



2J]. This full treatment is considerably more 



complicated and is presented elsewhere [30!]. For now we will only focus on the tensor modes, 
and do not consider the mixing of degrees of freedom. 

The main motivation for this work is to show explicitly that new effects can occur when one 
considers loop effects in a more general background then de Sitter space. We find new effects 
ffist of all in quasi de Sitter space, where due to the presence of an ultraviolet divergence one 
generates small, but physical corrections to the quantum Friedmann equations that cannot 
be subtracted by local counterterms. 

In section [IT] we briefly review our background geometry. In section IIIII we generalize 
the work of Ref. 20| and construct the massless minimally coupled scalar and graviton 
propagator in any e=constant space. In section |V] we calculate the one-loop effective action 
contribution to the quantum corrected Friedmann equations and renormalize the theory and 
in section |Vl] we study the associated dynamics in quasi de Sitter spaces. We conclude in 
section IVIIII We work in units: h= 1 = c. 



II. GEOMETRY 

The geometry we work in is the Friedmann-Lemaitre- Robertson- Walker geometry in con- 
formal coordinates 

gf,u = a^Vt^u ; Viiu = diag(-l, 1, 1, 1), (2) 

with the additional constraint, 

H TT ^ 

e = — T7^ = constant ; H = —. (3) 

H^ a 

Here a dot indicates a derivative with respect to cosmological time t, related to the conformal 
time 7] hj dt = adrj. The FLRW geometry obeys the Friedmann equations 

^-ip = ; -—-hp + p)=o^ K = 16nGM, (4) 

K Z K Z 

with Gn being the Newton constant, p and p are the energy density and pressure of the 
cosmological fluid. If one writes 

P = wp, (5) 



one immediately finds that ([3]) implies that w is constant. One can solve (jlj) for a to find 

(6) 



a{r])={ie-l)Hovf'^^' ^^ ; e = ^{1 + w) 






1-eH 
Notice that if e < 1, r^ is negative and if e > 1, 77 is positive. Hq is chosen such that the e — > 



expansion of H corresponds to the one given in 20| . An important geometrical quantity is 



y^y^^ = ^^i±^ = Xi-i\v-v\-^er + \\x-i\\'). (7) 

rjTj rifj 

Here the infinitesimal e > refers to the Feynman (time-ordered) pole prescription. In de 
Sitter space y is related to the geodesic distance I as y = 4sin^(|if/). If y < 0, points x 
are timelike related to x, and if y > 0, they are spacelike related. We define the antipodal 
point X of X by the map t] -^ —r]. Notice that, since in our coordinates rj is either positive 
or negative, this point is not covered by our coordinates. li y = 4, x lies on the lightcone of 
an unobservable image charge at the antipodal point x, see figure [TJ 







FIG. 1: The causal structure in the conformal coordinates ([2]). The plot assumes e < 1, so the 
coordinates (of an expanding universe) cover only the region 77 < 0. The wavy line at r/ = 
indicates future infinity. The lightcone of the point x is given by y = 0. If y = 4, the point x lies 
on the light cone of an unobservable image charge at the antipodal point x. 

The curvature tensors are given by 



^uVt^fs - S'^rif,^ 



7'\2 






where D denotes the number of space-time dimensions and a' = da/drj. 



III. SCALAR PROPAGATOR 



The construction of the graviton propagator in the geometry under consideration is very 
similar to the construction in quasi de Sitter space, as given in 20]. Therefore we will only 



give the main steps here. 

Since the graviton propagator can be expressed in terms of massless scalar propagators, 
we first consider the following Klein- Gordon equation for a massless scalar in D dimensions 



^^{U-iR)i^{x-x) = 

with some constant ^, where 
1 



-9 



D-^{D-l){D-2e)H^ 



iA(x: x) = iS (x — X 



n 



--d,. 



~gg''d, = -Ad'-{D-2)^d,^ , 



V-9 a 

denotes the scalar D'Alembertian. 

We now make the following Ansatz for the propagator 



(92 = ri'^'^d^d,] 



iA{x;x) = iahY^'-''/^\E{y). 



After the rescaling (TTTI) . the nonsingular part of the Klein Gordon equation reads 



,(4-,)(Af-(B(.-2))-|- 



- fl - e) 



' ((D -1){D- 2e)^ -^-{D-l){D-2)e+^{D- 2^ 



iE{y) = 0. 



This hypergeometric equation has a general solution 

^/ X . r. fD-l D-1 D y 

i^{y) = A 2F1 [-^— + i^D, — jyo; — ; ^ 

^ ^ /D-l D-1 D y 

+ B 2-ri h I/n, Vr)\ — ; 1 

^ ^V 2 2 2 4 



where 



V 



D 



D - 1\2 



:i-e) 



\d -1){D- 2e)i - \{D -1){D- 2)e + ^{D - 2)^ 



(9) 



(10) 



;ii) 



(12) 



(13) 



(14) 



The constants A and B are fixed by the singularity conditions. The vanishing of the antipodal 



34 



35 



36 



371], fixes A = 0. The constant B 



singularity (at y = 4), which leads to a- vacua 33|, 
is fixed by requiring that the Hadamard singularity at ?/ = sources the ^-function correctly. 
Notice however that there are values for vd where this is not possible. In particular if vd 
is half integer, larger or equal than 1/2, the hypergeometric equation is no longer a valid 
solution. For these particular cases one finds that one cannot source the 5 function correctly 
and remove the a-vacua. However, our solution is valid arbitrary close to these points. A 



well known examp 



in de Sitter space [38 



e o: 



such behavior is the massless minimally coupled (MMC) scalar field 



39|. 



The most singular term sourcing the 5-function is 

-( ^ r.fyy-''/' TiD/2)TiD/2-l) 

''^^^W^f^^-'^ (16) 

T{^ + un)n^-VD){\l~e\H,)''-' 

T{D/2) (47r)^/2 ■ ^ ' 

It is important to notice that, due to the rescaling ( TTTl) . B is indeed - as required - a constant, 

constituting a nontrivial consistency check of our Ansatz f ill I) . The MMC scalar propagator 

for a general, constant e reads 



usmg 



and ([6]), we find 



zA(x-x) =raa)-^(^/2 






2 

/D-1 D-1 D y 



Equation ( fTSl) is the generalization of the scalar Chernikov-Tagirov propagator 40|] to FLRW 
spaces with general, but constant, equation of state parameter w, and thus from Eq. ([6]) a 
constant e. 

IV. GRAVITON PROPAGATOR 

Next we consider the graviton propagator. The only difference from the analysis as given 



in 20| is in the coefficients 1^0,71- We consider the following action for gravity plus an arbitrary 



scalar field = 0(x). 

3 = 11 t/^xv^(i? -{D- 2)a) + J rf^xv^( - ^djd04>g"^ - V(0)) 

= 'S'eh + 'S'm 
and we write the metric tensor ^ as a background contribution g plus a perturbation h 

g^lu = 9t,v + V- (20) 



For the scalar fluid 



■p = (pit), the background energy density and pressure are 
PM = -lid<pf + Vi<P) = ^<p' + V{<P) 
PM = -lid<pf-V{<P) = ^<p'-V{<P), 



(21) 



where the background scalar field (j) = (j){r]) is a function of (conformal) time only. The 
corresponding tree level Friedmann equations are 



H' 



1 



-A- 



D-l {D-l){D-2) 



h+^^h^-Ia^ 



pM = 

Pm = 



2 2{D-2] 

dV 
b" + {D- 2)aH(j)' + a^— -(0) = . 



(22) 



Expanding the action up to quadratic order in h^j,,, and using the background equations 
of motion for gives ll|] 2G | 



^v 



^D 



d'^XXafSV^^V 



afj, Bu 



d' + ^Ha- + ^{D - 2)H-a-) {W. - \,X'' 



D~2 



4 



H' + H)aX5% 



(23) 



Xpo" ~r '-'o 



where we defined, 



''flU 



y Kd Ipfii/ y KCL X.^U-1 X.IJ.U X/1!^ nQliyX.-! X 9 Xf^l' 1 VPtlV ^ Xt^ly) 1 

(24) 



we added a gauge fixing term 



a 

2 



D+4 



V^(ai-^/2x^")lv,(ai-^/2~. 



(25) 



and where 5*08 denotes the second-order terms that vanish on-shell (c/. Eqs. ( 122|) ) 



On 



d^'xa'x^f, ( ^r^^V^ - ^V'-^V"'] 2{D-2)(H , ^ 



D-l 



-H'-- 



A kpm 



Xp<T 



2 2{D-2) 
+ j d^xa'x^p (-5i''v^^^%^ + l5lr5^,\''-^ {{D-2)H + ^{pM + Pm))xp. ■ (26) 

Even though these terms do not contribute to the graviton propagator, they do contribute 
to the one loop effective action and hence we must keep them. 
From the action 023p we find the graviton propagator for Xiiu 



^i-D/2~i-D/2^[^^^«/5] = (To),,°^zAo + (Ti),,"^Ai + (To),,"'3zA 



(27) 



where 



(To),."^ = 25("5f)- 



D-3 



v^"^ m) 



a/3 



per 



45?P^l?o'^ 



(T2 



0/3 



po- 



(D-2)(D-3) 
denote the relevant graviton tensor structures and 



(r^,. + iD- 2)8l5'^{t' + {D- 2)5m 



(28) 



(29) 



is the spatial part of the metric tensor. The prefactor a ' a ' in Eq. fl27j) comes from 
the fact that the x^iv held differs by a rescaling, x^lv = ci^^^^^'^'^^iv, from the V'^ti' fi^ld for 
which the propagator is calculated in |20| • The scalar propagators A„ for the pseudograviton 
i^fiu obey 

v^(n-n(D-n-l)(l-e)i72)«A„ = 2(5^(x-x'), (n = 0,l,2), (30) 



and thus they are given by flTSl) with 



V 



{D-lf , \{D - 1){D - 2)e - n{l - e){D - n - I) - ^{D - 2)e^ 



D,n 



+ 



4 ■ (1-e)' 

The ghost action associated with our gauge fixing is given by 



(31) 



5'ghost = / d^xri^pU'^ 



d' + ^Ha^ + jiD- 2)H\') 6-^ - {D - 2)a\H^ - H)5y, 



(32) 



where U is the ghost field, related to the ghost field of Ref. |2y] by a scale transformation 
U = aPI'^^^V . The ghost propagator is found to be 



jl^D/2~l-D/2^[^^p]^^. ~) ^ llP^k^{x-x)\%blblK{x-~x). 



(33) 



where the A„ propagators satisfy 



^ m - n[D - n - 1)(1 + e)if2 hA„ = ib"{x - x') , (n = 0, 1) 



Therefore the corresponding propagators are given by 0181) with 45 1 

.2 _[B- 1)2 \{D - \){D - 2)e - n{\ + t){D -n-l)-^{D- 2)e^ 



'D,n 



4 



+ 



'1-el 



(34) 



(35) 



Note that both the graviton and ghost propagators given above differ by a scaling 
a^^^/^a^^^/^ from the ones given in Ref. [20]. 



V. ONE-LOOP EFFECTIVE ACTION 



The one-loop effective action is defined as |4l| 



r = — 2(out, 0|0, in). 



(36) 



While in flat space and in the absence of external sources such a vacuum-to- vacuum transition 
can be normalized to unity, this is not possible in general curved space-times. Since our 
lagrangian is quadratic in the graviton and ghost fields, we can integrate them out to get 
the one loop effective action 



exp[«r] 



Vx^.u'DUVU exp 



exp 



'('S'eh + Sm) 



*(5'eh + ^M + 5*^ + 5*08 + Sghost) 

(•S'eh + 'S'm + 5*^ + 5'os + 5'ghost) 
det(^%) 



(37) 



where the step from the first to the second line can be made by noticing that the Jacobian of 
the transformation contributes as a D-dimensional delta function evaluated at zero, S^{0). 
In dimensional regularization such a term does not contribute. In fl37j) V and JF are the 
kinetic operators from the rescaled graviton (!23|l . including the off-shell contribution (126|) . 
and the ghost fl32l) . respectively, given by 



T^'^. = {d' + \{D - 2){D - 2e)H'a') {^6^^6: - \v,.v''') - (D - 2)(1 - e)H'a'6f/j;6, 



6D 



per 



flU 



-«Mo^W-7^'%02(Z}-2)(if 



D-1 



-H' - ^ + 



f^PM 



2 2{D-2 



+ a' 26lf6':X) + ^o^Vv,u ]{{D-2)H+ -{pm + Pm) 



K 



T\ = (d' + i(D -2){D- 2e)H'a')6-^ - (D - 2)(1 + e)H'aXS^ . 
From Eq. (1371) we obtain 

r = ^EH + ^M + '^Tr HVP\, + 5Vf\,] - iTi ln[^%] 
= -S'eh + "S'm + ^il ■ 



(3^ 



(39) 



While in principle one could - at least formally - evaluate the effective action, the object 
one is eventually interested in is the effective Friedmann equation, i.e. the equations of 
motion of the metric. Moreover in the present case there is the technical complication that 



10 



we need to work under the constraint that e is constant. As long as e remains small, there 
is no problem with imposing such a constraint in the equations of motion. On the other 
hand, imposing such a constraint in the action might change the dynamics substantially. 
By taking the functional derivative with respect to the scale factor a = a{ri), we obtain the 
Einstein trace equation, that is the —(00) + 3{u) component of the Einstein equation. Since 
in a FLRW universe there are only two independent equations, the second equation can be 
obtained by imposing the Bianchi identity. Thus our first equation of motion is given by 



6r 6{SnE + S, 



M) 



5a{l) 



Va^ 



5a{l) 
24 



+ 



5V 



IL 



K \ 



6a{l) 

-A + -if ) + 3pM - Pm 



6r 



IL 



(40) 



Sa{l) 



where V = J d^ ^x denotes the volume of space and pm and pm are the pressure 



and energy density associated to the matter action ^m and they are defined by, T 



M 



{2/ ^^—g)6S^[/5g^^ = —g,_iuPM — o-^^^^iPM +Pm)- We first focus on the graviton contribu- 
tion to 6rii^/6a{l). We first write the effective graviton action in terms of a: 



Eg [a] 



+ 



Tr In < ?7o(m'7^)(p'7<t)o 
1 



J V a 



d^ + UD-2) i(D-4)^ + - 



5V 



pa pu 



.(41) 



Now, since the term within the logarithm is just the kinetic operator, upon variation we will 
generate the inverse of this object. This inverse is of course the propagator. Taking the trace 
implies here both tracing over the indices, and evaluation at coincidence 4l|]. After taking 
the functional derivative we get, 

1 5V,[a] _ 1 5V'^[a] ^ 1 ^E^'H 



V 5a{l) V 5a{l) V 5a{l) 

where the latter term originates from variation of SVpcrfj.u in Eq. (|4T1) . We have 

1 ST'^ia] 



(42) 



V 5a{l) 



(«-2) (-^W?;')2^^ + (iw-i 



X 



D-2 /a 



a' d 
a? dl 



4 
AadP 



VpuV 



pa 



(43) 



^m{"^/3)^ Z 



aP^pa 



(x; x) 



11 



and 



V 5a{l) 



{D-2)\^{^-5i^5:^-\ 



^t^uV 



per 



X 



(D-3) 



a 



<o ?:{pk^) 



aO ^p<^ 



+ [^IW^' + -JIK)^ 









a' d 



■^ I n Q + n 77 I + 770 



a? dl 



1 t/2 



a(i/^ 



^M("^/3)>^^ 



aP'-^pa 



(44) 



(x; x) 



where 



a/3 ^ per 



(x; x) denotes the graviton propagator ( 1271) evaluated at coincidence. Eq. (l28l 



implies the following contractions, 









^M"^/3).^o.K,-) 



(^o)a/3pcr 
(T'lja/Jpo- 
{T2jal3pa 
{Tl)af3pcr 

{T2)a(3pcr 
{T2)a(3pcr 



D-1, 

1, 

D-1, 

2{D-3) 

D-2 ' 

4 

D-2' 



(45) 



75 



where {Ti)aj3pa = (^j)^/?^ V-ypVSa {i = 0,1,2). Other contractions vanish. After substituting 
l5l) into (jHj) one obtains 



1 Kl^] 1 
V^ (5a(/) 8 



L>(D-l)(D-2) 



D 



+-{D-l)iD-2) 



a'^ a" a' d\ I d^ ' 
a? a? a? dl) adP 

n _ 2 ) f — - — - —-) - - — 



a-" 



a^ 



a? dV a dP 



iAo(x;x)a ^ 



+ 



/2 // 

a a 



-(D-2?i- — - — — 

4^ Ha3 a2 a^d/ 



1 / \ 1 rf2 

- 3D-10 - — 
4 V / a (i/2 



«A2(x; x)a 



D-2 



(46) 



and 



1 ST'^Ja] 
V Sa{l) 



{D-2) 



+ {D-l) 



D{D-l) 



II /2 / 7 

a a ad 



(^-3)(^-^ + ^^j+aA + a^ 



X iAnfx; x)a 



D-2 



1^ 
adP 



(47) 



P-l)( 



// /2 

a a 



a' d \ . kV((1)) 
+ -^^) +aA + a- ^ ' 



+ 



{D+l) 



a? dl 



D-2 



II /^ I J 
a a ad 



T + 



a? dl 



kVU) 1 d^ ^ 

aA + a ^—^ H 

D-2 a dp 



iAi{x; x)a^~'^ 
«A2(x;x)a^"^ 
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This can be further siniphfied by making use of the on-shell relation (l22l) 

D-2 ^ ' 

For the ghost field we follow exactly the same procedure to obtain 

/2 



V 5a{l) ^ '{ 2^ 



(D 



a 



a-" 



a 



a^ 



a' d\ Id? 
a? dlJ adP 



+ 



1 ^n\(^ ^" ^' ^\ 31 d^ 



«Ai(x; x)a 



D-2 



(48) 

zAo(x; x)a^~^ 

(49) 



Upon combining fH6l) and (H9|) we obtain, 
1 i^r' „Jal f 1 



(B-2)( 



a' a" a' (i 



a^ a? dV a dl 



+ - 



1 d^ 



zAo(x; x)a 



D-2 



+ j(/>-i; 



f2 fi 

a a 






2Ai(x; x)a 



D-2 



1 

r 



a 

^2 



a' d 



'-jD-m--^--,j^]--„^^ 



3 1 rf2 
2arf/2 



+ 



-,3 



a' rf \ 1 (3D - 10) 1 rf 



«Ai(x; x)a 

2 



D-2 



1 /a' 

4^ ~^Ha3 a2 a2d/; 4 {D-2) adP _ 



tA2{x; x)a^-^ } , (50) 



where we made use of zAq = ^Aq. 

The next step is to evaluate the propagator at coincidence. From this point on we need to 
constrain our calculation to the case where e is constant. In this case the coincidence limit 
of the propagators are given by 

+ un,n)n^ 



tAJx:x) 



—-oP ^2A„(x;x) 
drj 

d^ 



\D-2ttD-2^^ 2 ) ^^ 



^D, 



(47r)T r(i + I/z),n)r(i-Z/D,n) 

Ha{D - 2)(1 - e)a^-^A„(x; x) 



.D-2„ 



dTf 



zA„(x;x) = H^a\D-l){D-2){l-eYa'''\An{x]x 



(51) 
(52) 
(53) 



where the v parameters for the ghost and the graviton are given by fl35|l and fl3T]) . respec- 
tively. Notice that the propagator is first evaluated at coincidence and only then hit by the 
derivative. Using these results in Eqs. ( l50l) and ( 1471) and Eq. ( HHl) we obtain 

1 



1 K+,h[<A 
V 5a{l) 



{D-l){D-2){l-e)H'a 



2„D-1 



-{D - 1){D - 2){D - A)eiAo{x; x) 



1 



--{D-l){D-2){2-e)tAi{x;x) (54) 

1 

~2 
1 

~4 



2(D + 2)- 


-3(1^- 


-2)e_ 


tAi{x;x) 


^4(D-3)- 


-{3D- 


-10)e 


2A2(x; X 
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and 



V Sa{l) 



{D - 2)H'a 



,n-ifD{D-l){D-2)e 



+ {D-1) D{D-l)-{D^-2D + 2)e 



D- {D-l)e iAo{x;x) 
iAifx; x) 



+ 



2D{D-l) - {'iD'^-QD+A)t+{D-l){D-2)e^ zA2(x;x 



(55) 



We substitute (1511) in (1541) and (135]) . add the two contributions and expand around D = A 
to obtain the nonrenormahzed one loop effective action [46| 
1 5Va+„M e(198-241e + 63e2)/i^-4 



g+gh[ 



H^a^-W 6a{l) 



4vr2 



D-4 



+ 



1 



167r2 



M - 1810e + 23076^ - 791e^ + 54e^ 



-2e(198-241e + 63e^ 



ij 



In 






8e(36 - 40e + Qe^ 

1 VI - 14e + e2 



54e(l-e)(2-e 



1 Vl - 14e + e2 



n2+ 2(1-.) )+n2 



2(1 - e) 



where '^/'(-z) = {d/dz)r{z) denotes the digamma function, and we made use of 



H 



D ,. D-4 17-4 



^""-'H' il + {D- 4) In {-)\+0 ((D- 



-4)^ 



(56) 



(57) 



where /x is an arbitrary renormahzation scale. 



Since e = —H/H"^ is in general a dynamical quantity, in order to renormalize the theory 
properly, one needs to subtract all divergent terms in Eq. (1561) containing powers of e. In 
order to do this we shall make use of the following counter lagrangian. 






9V(0) 



(58) 
where the last term denotes the Gauss-Bonnet term in FLRW spaces. This can be related 
to the standard form of the Gauss-Bonnet term by noticing that, since FLRW spaces are 
conformally flat and thus have a vanishing Weyl tensor, R^ivp^R^'^''" can be expressed as a 
linear combination of R^ and R^jR^'^ , 

2 



"'/ii^-' 



r> Dfiupa 

^""^ {D~l){D-2] 



(i?2 _ 2{D - m^.R^") . 



(59) 
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In order to fully renormalize the effective action Tg+gh (156|) nongeometric counterterms are 
required. These terms appear as a consequence of including the terms in the effective action 
'S'os ( l26l) that vanish on shell. The counter lagrangian ( l58l) is not unique. Indeed we could 



have chosen different counter terms |30l]. Since, based on the available information, there is 
no unique way to fix the counterterms, the form (l58l) of the counter lagrangian suffices for 
the purpose of this work. 

Varying the individual terms in the counter lagrangian fl58l) results in: 

^^ [ d^'x^/^R' = a^'-^H^i- 432e(l - e)(2 - e) 
V oa J \ 

+ 36(^4 - 34e + 356^ - Se^) (D-A) + o({D-4:)^, e 
^j^j d''xK^^gg^^''{d,(t>){d,<t))R = a^-iiJ^(l44e2(l - e) 

- 24e(2 - 8e + 5e2)(D-4)') + o(AD-Af, e) 

-16e2(4-e)(Z}-4)) +0((D-4)2,e) 
l^j-jd'>xy/^[R' - 3R,,Rn = a^'-'H' (36{1 - e)^) {D - 4) + o(p-4)^ e) , 

(60) 



where we used ([H]) and in the last step we again used the background equations of motion (l22l) 
and the following on-shell identities, 

v/^0' = V2p - 2)eaH ; |^(0) = 2{D - 1 - e)eH^ + 0{e) . (61) 

The divergent part of Eq. (l56|l cancels when the coefficients a, (i = 0,1,2) in the counter 
lagrangian ( !58l) are 

11 /x^-4 f 13 /i^-^ . 5 /i^-4 . , , 

where the af (z = 0, 1, 2) indicates a possible finite part of a^. 03 remains a free (infinite) 
parameter. 



15 



The renorinalized effective action FiLren is then obtained from 

1 l^riL.renM 



H^a^-^V 6a{l) 



^ I (/5o + /3ie + P2e' + Pse' + Pa^ 



167r2 ^ 
-2e(198-241e + 63e2) 



- 8e(36 - 40e + 9e^ 



4j 






54e(l-e)(2-e) 



,,1 Vl - 14e + e2 \ /I VI - 14e + e^ 

^'2+ 2(1-6) J+^b" 2a-e) 



(63) 



where the coefficients of the terms multiplying ey(167r^) {i = 0, 1, 2, 3, 4) are given by 



/3o 
A 

P2 



-48 + 5767r^(D-4)a3 



2168 



17287r' 



3 
4352 _„ 2 

1961 2 

+ 647r2 



8a^ + (Z}-4)a3 



396 



ln(47r)-7^ + 21n(/i/i?o) 



108a^ + 12a{ - 4^ + 9(D-4)a3 



482 



In(47r)-7s + 21n(/i/i?o) 



108a^ - 36a{ + 4^ - 9(D-4)a3 



+ 126 



54, 



ln(47r)-7^ + 21n(/i/iyo) 

(64) 



where ^0 is the expansion rate at which the \n(H/Ho) term in Eq. fl63|) vanishes. The 
formula (l63l) is one of the central results of our work. Even though (3^ in Eq. (!64|) seems to 
be fully specified by the one loop calculation, this is in fact not the case. Indeed, one can 
show that, upon adding to the counter lagrangian the counterterm C[ = a^^^/^R^d'^V / dcf^), 
P4 will become a function of 04, and thus unspecified. A similar statement holds for Pq: in the 
absence of the Gauss-Bonnet counterterm Pq has a definite value {Pq = —48). Since currently 
there are no physical measurements that specify the value of the Gauss-Bonnet counterterm, 
we conclude that 03 - and hence also Po - is unspecified by the one loop calculation (see also 



Refs. 



30 



m 



Other terms in Eq. ( l63l) . in particular the logarithm and polygamma functions, cannot 
be altered by local counterterms, and hence these terms constitute the physical graviton 
one loop contributions. According to the analysis of Ref. [30|], when mode mixing is taken 
account of, the poles of the ghost propagators coincide with those of the graviton, such that 
in the full analysis the digamma functions in the last line of Eq. (l63l) are absent (the same 
holds for the non-renormalized result (l56l)). 
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VI. DYNAMICS IN QUASI DE SITTER SPACES 



Equations PUI) and (!63|) . together with the Bianchi identity, give the quantum modified 
Friedmann equations. However, because of the complexity of (1631) . we will expand the 
correction in the limit of small e (quasi de Sitter space) [47|. 
When expanded in powers of e Eq. (!63l) gives, 
1 5rii,,,„[a] H"" 



Po - 252 



+ 



TT 

(3i + 374 - 792 In ( -j- 



7927i5 



e + 0(e^e)}.(65) 



Inserting fl65l) into Eq. fHOl) we obtain the following approximate Friedmann trace equation. 



H^-^+'-H 



/go-252 ^ 33 



ln(£ 



7e 



A + 374 



792 



{Pm - ^Pm) ■ (66) 



The quantum correction to the trace of the Einstein equation is the correction to the expec- 
tation value of the trace of the (quantum) Einstein tensor. 



6G = {n\5G\n) = -^^—^{n\5R\n). 



(67) 



From the symmetry of the background FLRW space we know that 6Gf^i, contains two in- 
dependent components: the first is the trace, and the second can be inferred from the 
corresponding Bianchi identity for 5G^v^ which is a consequence of the Bianchi identity for 
the background space Einstein tensor and of the covariant conservation of the matter stress 
energy tensor. Equivalently, one can view 5G^i, the 'stress energy' tensor {T^y)Q correspond- 
ing to the quantum corrections to fl66|) : then the symmetries of the FLRW determine its 
form to be, 

{T^u)q = diag(pQ, -pq, -pq, -pg) . (68) 

The covariant conservation of ( !68l) implies the following perfect fiuid-like conservation law, 

d 



dV ^^^ 
To solve for pq, we use the following Ansatz: 



3pq). 



/ H 



Xm + vH'H + {am + tH'H) In ^ + 0{t\ e 



which implies for the fiuid equation (169|1 that 
d 



^Pq) = a^ 4A/7' + {Av + 4A + a)H^H + {AaH" + (4r + Aa)H-'H) In 
at L V Hq 



H 



(69) 



(70) 



(71) 
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We read off (pg — Spq) from (l66l) and find that 



A 



647r2 

and thus 

1 r 

1 



(/3o - 252) , i; 



647r' 



:(/3o + A + 122 + 7927^), a = 0, r 



99 



87r2 



(72) 



(/5o-252)i7^ + /5o + A + 122 + 792^e)H'H - 792H'H\n 



f H 



\Ho 



+ 0{e',e) 



Pq+Pq = 



iPo - 252)H^ + ( l/?o _ /?, _ 458 - 7927^) H'^H + 792H^H In (J^)] + ^(e^, e) 



487r= 



■ (/5o- 252) /J2ij + 0(^2^ e). 



(73) 



The quantum corrected Friedmann equations become (c/. Eqs. f l22l) ) 



/5o - 252 



/5o + A + 122 + 7927s 



99 

TT 



F 

^ 



3H'-A+ ^" givff^ - ^""^^" " ^^ Gjvg^g + -G^i/^i/ln ( ^ 

Stt 87r 

+ 0{e\e)=8nGNPM 
Po - 252^ ^ 2 7V , xo^,2 



(74) 



2i/ 



Git 



-GnH'H + 0(e^ e) = 87rG,v(pM + Pm) • 



We shall assume that the matter contribution obeys an equation of state pu = wpM, with 
w constant. In this case we can combine the two equations as 



3H' + 



2H 



l + w 
where we defined 

A 

B 



A + AGnH^ + 



S H In ^ 



GNH^H + 0{t\k) 



0, 



(75) 



j3o-252 



in 



ilT 



Po - 252 - iPo + Pi + 122 + 7927^;) 



(76) 



-3(1 + w) 

Notice that, since we are free to choose (3o and /3i by a suitable choice of the coefficients Oq 
and as in the counterterms (see Eq. (1641) ) (the Gauss-Bonnet terms must be also included). 
Indeed, choosing /3o = 252 and Pi = —374 — 792je results in A = = B. In fact, A and 
B are not completely independent for a general value of w since from (1761) it follows that 
5 = 4A/[3(l + w)] + const. 

One can integrate Eq. fl74l) . The result can be expressed in terms of the roots of the 
quartic equation. 



A 



GnH"" + H^ 



A 



0. 



(77) 



In the case when A < all four roots ±H± are real, 



Hi 



2AGn 



AA 
-1±a/1 + — GjvA 



(7J 



The positive root H^ > corresponds to the one-loop corrected de Sitter attractor. From 
equation (1751) it follows that H+ is approached exponentially fast. More precisely, the late 
time limit can be approximated by the form. 



H = H. 



1 + 2 exp r - lit + 6q\ 



(79) 



where to order GmK and at late times Vtt^l, 



n 



m + w){i 



GnA 



A 
2" 



w 



B 



99 , / A 
— m — -TT 
27r \?,m 



(80) 



and 5q/VL represents an order G^A shift in time, which is unphysical since it can be absorbed 
in the definition of time. This means that quantum effects during quasi de Sitter phase induce 
an order Gj^K shift of the late time de Sitter attractor (which can be read off from H^ in 
Eq. 0781) ). At late times this de Sitter attractor is approached exponentially fast, with the 
characteristic time scale given by fi^^, which equals the classical time scale plus an order 
Gat A correction, as expected. In addition, there is an order Gat A shift 5q, which implies 
a time delay of 5q/VL. Note that 5q can be both positive and negative, depending on the 
sign and magnitude of A and B defined in Eq. ( !76l) . (The sign of 5q depends also on Hq, 
but a change in Hq can always be absorbed in a change m B.) This agrees with figure [21 
where we show if as a function of time both when 5q is positive (left panel) and when it is 
negative (right panel) (in the plots we have chosen A = and i? = 0). A positive (negative) 
correction 5q implies a greater (smaller) expansion rate H , and therefore a universe that has 
expanded more (less) before entering the late time de Sitter phase. 

At early times the quantum solution deviates more and more the classical solution, which 
approaches the Big Bang singularity at t = 0. Formally, the quantum one loop solution is 
not singular, and at large and 'negative' times (any negative time can be of course trans- 
formed to a p ositive time by an appropriate time shift) the solution approaches the quantum 



attractor 



42 



43 



44( 1 H -^ if_ defined in Eq. flTHl) . At this point the expansion rate becomes 



of the order the Planck scale, implying large higher loop corrections, such that this behavior 
cannot be trusted. 
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FIG. 2: Numerical solution to (|75p for H as a function of time. The red, dashed curve represents 
the classical behavior and the blue, solid curve includes our one-loop corrections. The late time de 
Sitter limit is clearly obtained. The quantum corrections lower the effective cosmological constant. 
In all plots we choose A = 3. This implies that we are effectively plotting the dimensionless 
variables: h = y/S/AH; r = y/A/3t and g = (A/3)GAr (A = 3 ; w = 1/3 ; Gn = 0.001 ; A = ; 
B = ; H{t = 0) = 10 ; Ho = 0.1 (left panel) Hq = 10 (right panel)). 

VII. DISCUSSION 



Before specializing the discussion to the two cases discussed above, we make some general 
remarks on the validity of our results. First of all, the correction we calculate is only valid 
when e is strictly constant. A nonconstant e would induce corrections to the propagators as 
calculated in section IIIII and unfortunately it is not yet known how to calculate these. This 
of course does not prevent one from using these propagators to calculate quantum correc- 
tions. One can then reasonably assume that, as long as in the final answer the change in 
e is sufficiently small, the error one is making is small and thus the results can be trusted. 
From figure [3] it is clear however, that there are regimes where e is far from constant and 
one should be careful to trust our results there. 

A second general concern is the issue of gauge invariance (invariance under infinitesimal 
coordinate transformations). When both a gravitational field and a matter field are present, 
the fluctuations in those fields are coupled and do not transform independently. Therefore 
one cannot self consistently quantize the gravitational fluctuations, without quantizing the 
matter fluctuations. The types of structures (e.g. the poles in the digamma functions), how- 
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FIG. 3: e as a function of time. The red, dashed curve represents the classical behavior and the 
blue, solid curve includes our one-loop corrections. The strong dependence of the behavior at early 
times on ^o is clearly visible. The parameters are: A = 3 ; tf = 1/3 ; Gn = 0.001 ; A = ; B = 
■ H{t = 0) = 10 ; Ho = 0.1 (left panel) Hq = 10 (right panel), see also figure El 



ever, are generic since they are naturally generated by the coincident limit of any propagator 
of the form ( ITSl) . and they do not disappear when matter fluctuations are included [30i] . 

Because the poles of the digamma functions ep = 1/2, 2/3, 3/4, ..,4/3, 3/2 that yield a 
divergent one loop effective potential (156|) are sufficiently distant from the quasi-de Sitter 
limit e ^ considered here, the results of our dynamical analyisis can be trusted as long as 
e -C 1/2. When this condition is satisfied, quantum effects do not change the fact that at late 
times the Universe asymptotes a de Sitter attractor, albeit with a modified expansion rate 
given by H+ in Eq. ( ITSll . The leading order quantum effect at late times has a contribution 
proportional to H^ to the effective energy momentum tensor. A contribution of this form has 

nnn 

also been found in earlier studies of graviton one- loop effects in de Sitter space [10| [25j [26j . 
The exact contribution is unknown because of the ambiguity in the counterterms. Depending 
on the choice of counterterms, the contribution could slightly increase or decrease the effec- 
tive late time cosmological constant. Although in our more general treatment, divergencies 
appear in the effective action, leading to the logarithmic correction to ( 1751) . these corrections 
have no significant effect at late times. 

At early times the contribution of quantum effects becomes more significant. However in 
this regime we have lost predictability, since the results strongly depend on the unknown 
part of the counterterms and the renormalization scale fi (Hq). Moreover, the assumption 
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that e < 1/2 and nearly constant appears to be violated. 

VIII. CONCLUSION AND OUTLOOK 

In this paper we calculated the quantum corrected Friedmann equations due to the one 
loop vacuum bubble from gravitons in a FLRW universe with constant e = —H/H^. The 
result has a divergence that contains terms proportional to eiJ^, t^H^ and e^H^, which can 
be renormalized using local counterterms, which include both geometric and scalar field 
counterterms. This is consistent with the result that in de Sitter space (e = 0) one loop 
effects lead to a finite constant shift of the cosmological term oc if^. We study the dynamics 
in the quasi de Sitter limit and find that they are not much different from the dynamics 
in true de Sitter space. Indeed, the quantum effects induce a shift in the effective, late 
time, cosmological constant ~ {A/3Q)GnH^, where A is an unknown parameter, that can 
be expressed in terms of the Gauss-Bonnet counterterm with an 0{1/{D — 4)) coefficient. 

Although our results are correct within the approximations used, the results described 
above should not be taken too literally. The propagators we used (and hence the singularity 
structure we find) are strictly speaking only valid when e = constant. Our analysis is correct 
as long as any time variation in e is small enough, which is indeed the case sufficiently close 
to de Sitter space. Indeed, our late time solution does have e — > 0. Therefore we have good 
reasons to believe that our solution approximates well the solution of the full theory, at least 
at late times and sufficiently close to de Sitter space. 

Another issue is that we choose our propagator such to describe a physically meaningful 
vacuum state. However, due to the evolution and mixing of modes, close to the de Sitter 
attractor the Universe will not be in a vacuum state, but in some excited state (that can 
be described by mode mixing in momentum space), which might infiuence our results. We 
postpone a study of this question for future publication. 

The next issue is the question of gauge invariance. Since there is both matter and gravity 
in our model, one should self-consistently take both fiuctuations in matter and gravitons into 
account. This issue is complicated due to the mixing of the degrees of freedom, and hence 
it is addressed in a separate publication 30|. Taking this mixing into account changes our 
results quantitatively, but since the singularity structure is inherent in the propagators (and 
those do not change), the logarithmic terms do not cancel, such that qualitative features of 
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the analysis presented here remain unchanged. 

Finally, an important question is what are the dynamics near the poles of the digamma 
functions e^ = 1/2,2/3,3/4, ..,4/3,3/2 where the one loop effective potential ( l63l) diverges. 
The corresponding dynamical analysis is performed in Ref. [30|]. Here we just note that the 
Universe typically gets stuck near the poles, such that each pole acts as a late time attractor. 
Probably the most important attractors are the two highest poles e^ = 3/2 and ep = 4/3 
(the latter is also the value of e in matter era). The latter pole is the late time attractor 



of a universe filled mostly with radiation and a cosmological term 30|, which represents a 



realistic composition of the Universe immediately after the Big Bang. 
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[46] Had we not included the contribution ()55[) from the off shell terms (126p in the effective action, 

instead of Eq. (f56l) ) we would get: 



5V 



g+ghV 



21/x^-4e(l-e)(2-e) 1-e 



ma^-W 5a{l) 



+6e(2 



7 In 



(1 



47r2(D-4) 
^2. ,H 



+ 



167r2 



Att 



14 In 



^ 



llE + 2^ 



1 



12 + Tie - 356^ -F 30e^ 



+ 2^- 



^, ,,1 , Vl-14e + e2 ^ , ^(\ VI - 146 + e^ 
''^'2+ 2(1-6) )+^(2- 2(1-6) 



1-6. 

+ 0{D - 4) . 



Note that the structure of the divergent term in this expression is simpler such that - unlike 
the divergence in Eq. (j56p - it can be renormalized by making use of the -R^ counterterm alone. 
[47] The question of the dynamics around those e for which the one-loop potential ()63p exhibits 
poles is addressed in the companion paper [3C|. 



